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Using 'perturbation theory, an expression jor the variance of the power 
of each mode of a multimode waveguide with randomly-coupled modes is 
derived. The variance builds up from zero to a constant value as a function 
of z {length along the waveguide). For most cases of interest, the variance is 
equal to the square of the average power. This means that the power of each 
mode of a system of randomly-coupled modes of a multimode waveguide 
fluctuates like tfie short-term time averaged power of a narrowband electrical 
signal the voltage of which is a random variable with Gaussian probability 
distribution. 

I. INTRODUCTION 

The behaviour of waves propagating in multimode waveguides can 
be described by coupled equations for the amplitudes of each mode. 1 
This description is rigorous, but has the disadvantage that the coupled 
wave equations usually cannot be solved. It has been shown that a 
much simpler description is possible if we limit our interest to knowledge 
about the average power carried by each mode. 2-4 Coupled equations 
for the average mode power have been derived and applied to the 
problem of wave propagation in multimode dielectric waveguides. 4,5 
However, the description of multimode waveguides in terms of average 
power is incomplete unless some information is available about the 
fluctuations of the actual power about the average value. With the help 
of the same perturbation approach that was used to derive the coupled 
power equations, 4 we derive in this paper a differential equation for the 
variance of the power. 

The result of our perturbation theory is expressed in terms of the 
cross-correlation and the average power of the modes. In order to 
evaluate this expression, we need to make several assumptions. It has 
been shown in an earlier paper 5 that the average power settles down 
to a steady-state distribution of power versus mode number that is 
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independent of the initial excitation of the waveguide. Using this concept 
of the steady-state distribution and the further assumption that the 
cross-correlation between the modes is small, we can solve the differential 
equation for the variance. We find that the variance builds up from zero 
values at z = 0, to a constant value which is equal to the square of the 
average of the mode power. The relative fluctuation of the power of 
each mode is thus 100 percent. This means that the power in each of the 
randomly-coupled modes behaves like the short-term time-averaged 
power of a narrowband electrical signal the voltage of which is a Gaussian 
random variable. 

II. DERIVATION OF THE DIFFERENTIAL EQUATION FOR THE VARIANCE 

Our starting point is the set of coupled wave equations for the slowly 
varying wave amplitudes (envelops) A which are defined by 

a, = A,e-^' (1) 

with a, being the rapidly oscillating mode amplitude. The coupled 
wave equations can be expressed in the form 

^= E^pV 1 *-*""-" (2) 

dz M =i 

z' is used as a convenient reference point. The parameters /3 V are the 
propagation constants of the modes. The coupling coefficient can be 
expressed as a product of a constant term times a function of z. 

c„ = KJ(z). (3) 

If we define 

K„ = 0, (4) 

we can drop the restriction y. j* v in (2). Conservation of power leads 
to the relation 4 (the asterisk indicates complex conjugation) 

K„ = -K* . (5) 

The perturbation theory uses the approximate solution of (2) 

AM = A,{z') + E KJLtf) f Kx)e iU) '->» )z dx. (6) 

The power of mode v is 

P,(z) = | a, | 2 = e-' ( -"' | A, | 2 . (7) 

a v = —2Im0 v is the power attenuation coefficient of mode v in the 
absence of coupling. Throughout our derivation we assume that the 
losses are so slight that we can approximate exp ( — a v (z — z')) by unity. 
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We use the loss term in (7) only to modify our equations for the lossy 
case. 

The variance of the power of mode v is defined as 

(AP,) 2 = (PR - PI , (8) 

with the simplified notation 

A = <P.>. (9) 

The derivative of the variance can be written with the help of (7) 
(replacing the exponential term by unity) 

I|(^ = {(^^^) + c.c.}-P.f. (XO) 

The expression c.c. indicates that the complex conjugate of the first 
term in the bracket must be added. The derivative of the average power 
has already been evaluated so that we do not need to express it in terms 
of the wave amplitudes. With the help of (2), (10) can be written as 
follows: 

|| (A/\) 2 = { £ {A,A*A,A*Kz))K,/"'-°> u + c.c.) - P, &- (11) 

We now follow the technique that was developed in Ref. 4. We replace 
all the amplitudes in (11) with the approximate solution (6), but keep 
only terms up to second order in K VII . The first-order terms vanish 
if we assume that j(z) is statistically independent of A(z'). This assump- 
tion is justified if we let z — z' be much larger than the correlation 
length of f(z). For the same reason, we write the ensemble average of 
products of the field amplitudes with terms containing /(z) as a product 
of an ensemble average containing only amplitude terms, times an 
ensemble average of a term that contains only f(z). We thus obtain 

- {Z [K„K, i (A*AtA>A!)e" 2 ''-''-' t) ' 
•/' (/(2)/(.x-)>e , ' ( "" / " ,(l " I) dx + 2K„K* i {AMA li A*)e im -^ u 
•/" (Kz)Kx))e- i<fi '- piUl - ) dx + K„K llt (A.A*A s A*)e ilfi '- f ' ,u 
'f </(«)/(aOy ( ''~ / ' l)(, ~' J efa] + c.c.} - P, £ h.Jp, - P,). (12) 
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The amplitudes A< are understood to have the argument z' . The last 
term was obtained by using the lossless coupled equations for the average 
power derived in Ref. 4. The power coupling coefficient is defined as 4 

K = I K » I 2 F (ft- - ft). ( 13 ) 

The power spectrum F% - 0„) is the ensemble average of the absolute 
square value of the Fourier coefficient of f(z). 

The next step in the derivation is based on the realization that only 
nonoscillatory terms contribute appreciably to the growth of the 
variance as a function of z. We thus neglect all but the nonoscillating 
terms in (12). This procedure is reinforced by the fact that the ensemble 
averages of cross terms of amplitudes are likely to be smaller than the 
ensemble averages of absolute squares of the amplitudes. The first 
term in (12) causes some concern since it appears that there may be 
several combinations of n and 8 in addition to n = 8 = v that contribute 
nonoscillatory terms. However, the uneven spacing of the modes along 
the axis makes it appear unlikely that combinations of modes can be 
found for which the exponent of the exponential function in front of 
the integral vanishes. However, even if a few such combinations could 
be found, we could still consider the term belonging to such combinations 
as small because of the lack of correlation between the amplitude coeffi- 
cients belonging to different modes. Since K^ = 0, we find that the 
first term in (12) does not contribute appreciably to the derivative of 
the variance and can be neglected. The integrals can be expressed in 
terms of the power spectrum of the function /(z) as was shown in Ref. 4. 
We thus obtain, with the help of (5) and (13), and dropping all oscillatory 
terms 

T (AP,) 2 = 2 £ KMP^) - P.P, - (Ai\) 2 ]. (14) 

az ^=i 

Assuming that P,(z') « P,(z), we use z as the argument of P, . Finally, 
we introduce losses into the theory. The example of the coupled power 
equations serves well to illustrate the procedure. Neglecting losses we 
obtain 4 

^r - 2 KJlP, - A). (is) 

az M =i 

According to our derivation, we have used the approximation P, = 
( | A„ | 2 ). Using (7), we obtain 

*4H*- (•*•+$ >""""• (16) 
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The left-hand sides of (15) and (16) are identical according to our 
derivation. By substituting (16) and assuming again that a(z — z') <<C 1, 
we obtain 

«,P, + ^ = Z Mft - ft). (17) 

Equation (17) is identical with equation (29) of Ref. 4. There, we 
introduced the loss simply as a phenomenological parameter. Our 
present treatment shows how the loss term can be obtained directly 
from the derivation based on perturbation theory. By applying the 
same reasoning to (14), we obtain 

| (AP,) 2 = -«,(AP,) 2 + 2 £ fc,„[2<(P, - ft)(P„ - ft)) + ftft]. (18) 

The parameter k, is defined as 

k k = 2a, + 2 £ fe,, . (19) 

We regrouped the terms under the summation sign in (18) in order to 
express (P,P») in terms of the cross correlation ((P„ — P v ){P ll — ft)). 
Integration of (18) yields, finally, the desired expression for the variance 
of the mode power 

(AP,) 2 = (AP,) a ,„„ + 2e-" 



/: 



e" 1 E *m[2<(P. - ft)(P, - ft)) + P,ft] rf.r. 
».=i 

(20) 

Equation (20) is the solution of the variance problem. The expression 
in brackets under the summation sign can be positive or negative, so 
that the variance can increase or decrease with increasing z. 

III. EVALUATION OF THE VARIANCE FOR SPECIAL CASES 

In order to be able to evaluate the general expression (20) for the 
variance, we would need to know the cross correlation and the average 
power as functions of z. The average power can be obtained by solving 
the coupled power equations. However, the cross correlation is not 
known. It appears reasonable to assume that the cross correlation may 
be small in many cases of practical interest. One would not expect to 
obtain small values of the cross correlation for only two modes because 
as one mode gains power the other must lose an equal amount of power. 
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However, for large numbers of modes, it appears reasonable to expect 
that the cross correlation between different modes may be small. 

It is known from the theory of coupled power equations that the 
distribution of power versus mode number settles down to a steady 
state. 5 Once the steady state is reached, each mode decays with the 
same attenuation coefficient. The shape of the distribution of average 
power versus mode number remains unchanged, but its level decreases 
exponentially with a power attenuation constant a, . If we launch a 
power distribution at z = that corresponds to the steady-state distri- 
bution, we obtain power averages that do not change with z except for 
a common exponential decay term. Assuming, therefore, that the cross 
correlation is negligible and that the steady-state power distribution 
is launched into the guide, allows us to solve (20) immediately. Using 
P, = P„e~ a " we obtain for (AP,tf. fl = 

g4 L -^-^fc} (21) 

with k, given by (19). Equation (21) represents the relative fluctuation 
of the power of mode v. It shows clearly that the relative fluctuations 
build up from zero to a constant value which is reached when the 
z-dependent exponential function in (21) becomes negligibly small. 
The shape of the steady-state power distribution depends on the interplay 
between the coupling between the guided modes and the loss of power 
to radiation. The loss coefficient a, that appears in the coupled power 
equations (17), depends on the mode number. Usually higher-order 
modes suffer more losses than lower-order modes. Modes with a large loss 
coefficient carry only little power once the steady-state power distri- 
bution is reached. Modes with small average power are of little interest. 
Concentrating on those modes that carry appreciable amounts of power 
allows us to neglect the attenuation coefficient a, that appears implicitly 
in (21) through relation (19). If all the guided modes couple strongly 
to each other, they are also strongly coupled to the radiation field, thus 
losing a large amount of power by radiation. Since reasonably low loss 
operation is of most interest in practical applications, we can limit our 
discussion to the situation where only neighboring guided modes are 
coupled to each other. This means that h y)l is small for large values 
of | n — v | . The sum over h vll thus extends only over those values of n 
which are close to v. The steady-state power distribution is continuous 
in the sense that neighboring modes carry nearly equal amounts of 
power. Neglecting the small steady-state loss coefficient a. compared 
to the sum over the coupling coefficients h fli , and using the fact that 
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the power of neighboring modes is nearly equal, allows us to obtain 
for z — > oo from (21) the important relation 

7F 5 « 1- (22) 

IV. DISCUSSION OF THE RESULT 

The relative fluctuation of the power of those modes that carry 
appreciable amounts of power is approximately 100 percent. Such 
fluctuations are not unusual, however. The short-term time-averaged 
power carried by a narrowband electrical signal, the voltage of which 
is a Gaussian random variable, is known to fluctuate in the same way. 
The probability distribution for P, can, in analogy to the electrical case, 
be assumed to be 



W(P r ) = ±exp(-^)- 



(23) 



From this analogy we can immediately state that the relative fluctua- 
tions of the power of M modes (assumed to be uncorrelated) is equal 
to M~ u \ 

The fluctuations that we are considering do not occur in time at the 
output of any given waveguide. They are fluctuations of random vari- 
ables in an ensemble sense. If we were to measure the power in a given 
mode for each of a large number of similar waveguides, we would 
expect to obtain results that fluctuate according to (22). Equation (22) 
thus tells us the accuracy of predicting the value of the power in a given 
mode on the basis of the coupled power equations. Since it is very hard 
to measure the power carried by one individual mode of a multimode 
waveguide, we are more likely to observe the power P M in a fairly large 
number of M modes simultaneously. In this case, we expect to obtain 
fluctuations according to the law 

t - W 

It is helpful to remember that the power of all N modes does not fluc- 
tuate at all. 

V. CONCLUSIONS 

We have discussed the problem of the relative fluctuations of the 
power in individual modes of a multimode waveguide in the case that 
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the modes are coupled by a random coupling function. Our discussion 
was limited to the c.w. case and does not directly apply to pulsed 
operation. We derived the general expression (20) for the variance of 
the power in terms of the cross correlation and the average power carried 
by the modes. Under the assumption that the modes are approximately 
uncorrelated among each other, and assuming further that only neigh- 
boring modes are coupled, we found that the relative fluctuations are 
nearly 100 percent. This result is reminiscent of the fluctuations of the 
short term time averaged power of a narrowband electrical signal the 
noise voltage of which is a Gaussian random variable. 

Cross correlation between the modes can either increase or decrease 
the fluctuations depending on the sign of the cross correlation term 
in (20). It is reasonable to assume that the sign would tend to be nega- 
tive. As stated earlier, coupling between only neighboring modes is 
necessary for low loss operation. If mode v should, at a given point on 
the z axis, carry more than the average amount of power, we conclude 
that this power has been transferred from the neighboring mode (or 
modes) n so that this mode is expected to have less than the average 
amount of power. The sign of the two factors in the cross correlation 
term must thus be different so that the term assumes a negative sign. 
This qualitative discussion indicates that correlations among the modes 
would tend to reduce the variance (AP„) 2 . I have observed fluctuations 
of the mode power as large as those predicted by this theory in numerical 
solutions of coupled line equations with random, band-limited coupling 
function. This "experimental" result confirms the assumption that cross 
correlation between modes does not appreciably reduce the variance 
of the power fluctuations. 
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